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Abstract 

This paper discusses two methods of developing models for the rapid pressure-strain correlation 
term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One 
is a perturbation about isotropic turbulence, the other is a perturbation about two-component 
turbulence — an extremely anisotropic turbulence. A model based on the latter method is proposed 
and is found to be very promising when compared with DNS data and other models. 

1. Introduction 

Transport equation models for second moments such as the Reynolds stresses are gradually 
being accepted as advanced tools for studying complex turbulent flows. This is due to the fact 
that important turbulent physics can be represented by the terms in the second moment equations, 
e.g., the pressure-strain correlation tensor, the dissipation rate tensor, etc. The correct modeling 
of these terms will enable us to capture important turbulent physics and to improve the prediction 
of complex turbulent flows. 

Over the past severed years, considerable attention 1-7 has been directed towards the develop- 
ment of models for the pressure-strain correlation term which is considered to be one of the most 
important terms in the second moment equations. Many attractive models have been developed 
by using a variety of theories and principles, such as invariant theory, rapid distortion theory and 
realizability, etc. The present paper will concentrate on the rapid part of the pressure-strain corre- 
lation term. We will first describe some of these models and related theories. Then we will explore 
whether or not the intended behavior of the models is possible by examining the results of direct 
numerical simulations (DNS) of homogeneous turbulence?. By analysing these models using DNS 
data, we are able to identify the deficiencies of the models and the sources of the deficiencies. We 
are also able to identify the applicable range of flows for each model and provide the direction for 
improving the model’s performance. For example, in the Launder, Reece and Rodi model 1 , we 
are able to suggest a function form for their model coefficient Ci which leads to good agreement 
with the most important shear component for all available DNS shear flows. In this paper, a rapid 
pressure-strain correlation model based on Shih and Lumley’s method 3 has been developed using 
a perturbation series about a two-component turbulence state. 


2. Modeling of the pressure-strain correlation 

The Reynolds stress equation for homogeneous turbulence in terms of the anisotropy tensor 
b XJ can be written as 


d h 

dt *' 3 
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where b X3 = v^uj/2k — 8 XJ /3 is the anisotropy tensor, U X>J is the mean velocity gradient, k and e are 
the turbulent kinetic energy and its dissipation rate. Ily and e X} are the pressure-strain correlation 
tensor and the dissipation rate tensor, respectively. The latter are the new unknown terms in the 
Reynolds stress equation and must be modeled. In this paper we only consider the rapid part of 
the pressure-stain correlation term 11^. Its exact expression is 

ng>-a P) 

where 

Xpjqi = ~ ^ J J J [u 9 (r)« t (r , )] >P j ^ dv (3) 


2.1 Phenomenological modeling of the rapid pressure-strain correlation 

Equation (1) indicates that the relationship between the pressure-strain correlation and the 
anisotropy tensor b xo should in general be nonlinear due the nonlinear term 2b X] b pq U p ^ q . Equation 
(2) reveals that modeling of the rapid pressure-strain correlation requires a model of the fourth 
order rank tensor X pjqt . In the literature, X piqx is assumed to be a function of b XJ . However, 
Eq. (1) indicates that X pjqx should at least be a functional of b XJ , i.e., it should depend on the 
history of b XJ in the whole flow field. Therefore, the real problem is quite complicated. Here, 
as an engineering approximation, we will assume, like most other researchers, that Xp ]qx is only 
a function of b x3 at local point in space and time. We notice that this assumption will at least 
lose the effect of dimensionality on X P3qx and may cause the problem in modeling of rapid rotating 
turbulence as pointed out by Reynolds 9 . However, in practice, this assumption is quite appropriate 
for many turbulent shear flows. Furthermore, we expect that X P3qi should at least be a quadratic 
function of b XJ because of the nonlinearity between 11,^ and b x3 . In fact, from a method of invariant 
theory of rational mechanics we may find that the most general tensorial function form for Xp 3qt 
in terms of b X3 is 3,6 

2^ ~ oc\8 qx 8 P3 d* c^2 -f- 8 q j8 p i) d- &3 8q%b p j d- U4 8 p j b qi 

d- U5 (fipqpxj d~ 8 X j b pq d- 8 3q b p 1 d- 8 p i b q3 ') d* OtfiSqibpj d- ^lb P3 b ql 

d- ot%{8 pq b 33 d- 8 X3 b pq + S 3q b pt + 8 px b q3 ) + ot§b ql b P3 
+ o‘io(bp q b %3 d- b q3 b px ) d- an6 gt 6pj + 0:12 b P3 b 2 qt 
d- Ui3 ippqb %3 d- b X3 b pq + b q3 b pi + b px b q3 ) 
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Equation (4) already satisfies the symmetry properties of Xp 3qx = X 3pqx and X P3qx = X p3iq required 
by Eq. (3). The number of coefficients in Eq. (4) can be reduced to nine by the following basic 
mathematical properties required also by Eq. (3): 


X ppqx — 2k(b qx d” X p3q3 — 0. 


(5) 


When the term ( X P3qx d- X pxqj )U Ptq is formed, the coefficients will be further reduced to seven as 
shown by Johansson and Hallback 7 because of the following tensorial identity relation^’ 7 : 


( b 2 xp 8 3q 


+ b] p 8 xq - b 


f zj bpq 


+ b t qbjp 


‘ bpqbij )Spq — 0 


( 6 . 1 ) 


2 



( 6 . 2 ) 


( b 2 p b jq + b%b tq - b tJ b 2 pq - b%b pq + \lII8 xp 6 3q )S pq = 0 
[ b tq b 2 p + b 2 q b 2 p — 2 b XJ b 2 q — II(b X3 b pq — b xq b 3P ) 

+ — III(b tp 8 3q + b-,p6 tq — b pq 8 %3 ) ]S P9 = 0 (6.3) 


2.2 Truncated anisotropy power series models 

The seven undetermined coefficients in Eq. (4) axe in general functions of the invariants of 
b t3 : II = —b tj bj t / 2 and III = b X] b 3 kbki/%- Because the absolute values of II and III are always 
less than unity, one may expand the coefficients in a power series of II and ill 6,7 : 

a, = C t (1) + C, (2) II + C t (3) III + C, (4) II 2 

+ C, (5) II III + Cf > III 2 + • • • (i = 1, 15) (7) 

The coefficients in Eq. (7) are independent of II and III hence are considered to be constant. 
Therefore, if their values are determined by using a particular flow, then the values should be valid 
for other flows as well. That is, Eq. (7) (with an infinite number of terms) will provide us with 
a universal rapid pressure-strain correlation model. However, in practice, we must truncate the 
series to a certain power of II and III and use certain flows and constraints (such as homogeneous 
flows and realizability) to determine the coefficients in Eq. (7). Having done that, we then must 
verify if such determined coefficients are constant and are valid for other flows as well. 

A first order power truncation model — LRR model 

If we decide to only keep the terms in Eq. (4) which are of up to order 0(b) (the norm of h 3 ) 
and note that II = 0(t?) and III = 0(b 3 ), then the coefficients a^, ■ • • , <*15 must be set to zero 
and in Eq. (7) we only need to keep the first term C^\i = 1,5) which are constant. In this case, 
Eq. (4) becomes 


2k 
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Using the properties in Eq. (5), there will be only one undetermined coefficient left and the 
corresponding rapid pressure-strain correlation model Eq. (2) becomes 
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where S XJ = | (U t>3 + U 3fX ) and Q tJ = |(U t)J — U 3 , x ). Equation (8) is the Launder, Reece and Rodi 
(LRR) model 1 . The undetermined coefficient Ch is taken to be constant and was set at a value of 
0.4. 

Now, we may use DNS data of turbulent shear flows 8 to examine the coefficient C 2 in Eq. (8). 
For homogeneous turbulent shear flows, there are only three independent components Hu, H22 
and Hi2- From Eq. (8), each component of H,j can be used for determining C 2 . Figure 1 shows the 
values of C 2 deduced from the DNS data of Hu, n22 and ni2 for the case of C128W. Apparently, 
C 2 is not a constant. Different component gives a different value of C2. This situation is also ture 
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for other flows (C128U, C128V, and C128X). The above result indicates that the first order power 
truncation model can not be expected to work for all the turbulent components. However, Figure 
1 shows that for the 1-2 component, C 2 = 0.4 is indeed a good approximation. This indicates that 
LRR model is probably good for two dimensional turbulent shear flows since the 1-2 component 
dominates the flow. However, if the flow is dominated by normal components, then we cannot in 
general expect LRR model to give a good prediction. 

In order to improve the model performance, a natural way is to pursue a higher order power 
truncation. 


A fourth order power truncation model — JH model 

Now if the terms in Eq. (4) axe retained up to order (6 4 ), then it will require the following 
truncated power series in Eq. (7): 

= C t (1) + C t (2) II + C t (3) III + Cf 4) II 2 (i = l,2) 
a, = C t (1) + C[ 2) II + C t (3) III (2 = 3, 4, 5) 
a t = C t (1) + C, (2) II (2 = 6, 7, 8, 9, 10) 
a, = C; (1) (2 = 11, 12, 13, 14, 15) 

In this truncation form, there are 32 ’’constants” d^\c[ 2 \ • • • , need to be determined. Using 
the constraints in Eq. (5) and (6), and the condition of realizability 3 : 

Up y qXp a q a = 0 if = 0, (9) 

Johansson and Hallback 7 found that there are only four undetermined constants left and the rapid 
pressure-strain correlation model becomes 


where 
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and 

Bi = fo - 8 foil + 24 frill, B 2 =fo- 8 foil 

Bz = fo — 8 foil, Bi = fo, B 5 ~ fo, Bq = foo 

and 

fo = -^(3 + 60/? x + 48 fo - 40 fo), fo = ~\~ 132/? 2 + |/?io 

fo = 2 + 60(/?2 + /?3)j fo = g/?l - ^(/?2 + /%) 

= — gg(g + 21/?2 + 10/?3 + — /?io)> fo = “ 220^8 + + ^ 3 ) 

The coefficients fo,fo,fo and foo are undetermined constants. To determine these constants, 
Johansson and Hallback 7 used rapid distortion theory to analyze an irrotational strain imposed 
on initially isotropic turbulence and a pure rotation imposed on anisotropic turbulence . They 
suggested that 

fo = fo — 0.0295, fo = -0.0484, fo 0 = 14 

This fourth order power truncation model Eq. (10) was shown to be in excellent agreement with 
irrotational RDT results. The result indicates that for the RDT irrotational flows, the fourth 
order power truncation is appropriate and the truncation error is quite small. This model is also 
in excellent agreement with DNS data of axisymmetric contraction flows shown in Figure 8. 

Now let us use DNS data of turbulent shear flows 8 to examine the “constants” of fo,fo,fo 
and foo. Since only three components of II tJ are independent, Eq. (10) can determine only three 
coefficients. We set fo = —1/7 and use the DNS data of Du, n 2 2 and II 12 to deduce the coefficients 
fo , fo and foo. Figures 2 shows the values oifo,fo and foo deduced from the DNS flows of C128U, 
C128V, C128W and C128X. Apparently, these coefficients are still not constant. If we believe that 
the DNS data are realistic and that the formulation of Eq. (10) is correct, then the above result 
indicates that the fourth order power truncation still have too large of a truncation error for the 
flows considered above. No constant values of these coefficients can ensure that the model will 
have good performance. To overcome this problem a much higher order truncation then must be 
retained to reduce the truncation error. However, at the present time, it is not clear what order of 
truncation is needed for a realistic turbulent flow. In addition, it is noticed that the values of fo,fo 
and foo deduced from RDT results 7 are very different from those deduced from DNS flows 8 . This 
inconsistency raises a question as to which of these flows are more realistic and can be considered 
as benchmark for turbulence model development? 

2.3 Truncated F power series model 

Here we discuss another truncation method which is a perturbation about the two-component 
(2-C) turbulence state. The deviation from the 2-C state can be measured by a parameter F which 
was introduced by Lumley 2 : 

F = 1 + 9 11 + 27 III (10) 

From the definition, F will vanish if turbulence is 2-C (in which the turbulence becomes extremely 
anisotropic) and F will equal unity if the turbulence is isotropic For a realizable turbulence, F is 
always positive and varies between zero and one. 

In modelin g the rapid pressure-strain correlation term, the focus should be directed towards sit- 
uations where the turbulence anisotropy is strong (F « 1) rather than cases where the anisotropy 
is weak. This is because the Reynolds stress transport equation mainly emphasizes the effect of 
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anisotropy. For a flow with strong shear, F could be driven to zero while b is about 1/3. Based on 
the above consideration, we propose to expand the coefficients in Eq. (4) in a power series of F: 


a t = cj 0) (III) + eg (III) F + eg (III) F 2 (i = 1, 15) (12) 

where the coefficients eg, eg, eg, • • • are functions of the third invariant III of b XJ and they can 
be expanded in a power series of III : 

eg (III) = eg + eg III + eg III 2 * + • • • (12a) 


For small III, as a first order approximation they can be considered as constant. We will examine 
the behavior of the corresponding rapid pressure-strain model. 

Let us now start with the most general form, Eq. (4). Using the conditions in Eq. (5) we may 
obtain six relations between the fifteen coefficients which will reduce the number of coefficients to 
nine. As a result, we obtain 


cc l 
C*2 
<*3 
a 4 
a 7 

OC 9 


2 II , n , III , , e , 

— + — (4<*6 — 2 a$) + —(an + a 4 2 — 6013 ) 
lo 0 0 

_ 30 — jf( a6 — 3c * 8 * * * ) — !o-( 3ai1 + 3c * 12 + ^ ai3 ^ 

1 11 II, n 0 , III , n . 

— - — —0:5 + -^-(0:11 + 3 ai 2 + 8013) ^-( 3 ai 4 + oqs) 

1 4 2 II, n , 2 

3 " 3 0:5 + T^ 11 + 2 “ 13 ^ ” 3 II/ai5 

2 4 211, 

~ jj<*10 — g<28 + -g-(Q!l4 + ai5) 

11 1 II, _ . 

— 3 " as ~ gttio - <*6 + y (<*14 + 7ai5 ' 


(13) 


Equation (13) will ensure the rapid pressure-strain correlation model to satisfy the basic mathe- 
matical conditions in Eq. (5) while the values of the nine coefficients: 05 , ae, as, quo, an, 0 : 12 , 
< 213 , <*14 and CK 15 can be arbitrary. In addition, the realizability condition of Eq. (9) will provide 
three additional relations for the c/s at the 2-C state, i.e., F = 0. They lead to 


<*5 = - Jq + ^(<*11 + ai2 + 3 4 * <*i3 - 10aw/9) 

+ y ( 3 aii + 3 ai 2 + 8013 — 2a 15 — 60:6) 

1 1 

<*8 = g<*13 - g<*15 

3 1 

aio = + — (9an + 9ai 2 + 66a 13 - 10ai 5 ) 

+ Y^(9an + 9ai2 -I- 24ai3 - 6a 15 - 18ae) (14) 


Equation (14) is valid only at the 2-C turbulence state and will ensure the rapid pressure-strain 

correlation model satisfies the necessary realizability condition of Eq. (9). Therefore, we conclude 

that Eqs (13) and (14) will ensure the model satisfies the conditions in both Eq. (5) and Eq. 

(9) with arbitrary values of the six coefficients: 05 , 0 : 11 , 012 , <* 13 , Oi 4 and a 15 This will allow us 

to set any of these coefficients to zero in order to simplify the model while maintaining its basic 

properties described by Eqs (5) and (9) 
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Tensorially quadratic model 

Here we will set all the six undetermined coefficients a$, • • • , ais to zero to obtain the most 
simplified model. The necessity of retaining these coefficients will be left for a future study. Under 
this consideration, Eq. (13) becomes 

2 2 II 1 2 II 1 11 

a e , a2 = -- + — a 8 , as = -j- T a 5 , 

1 4 2 4 11 1 

0-4 = - - -<* 5 , (*7 = --aio - -as? <*9 = — - g^lO- (15) 

and Eq. (14) becomes 

<*8 =0, a-io = — ^ if F = 0 (16) 

For F 7^ 0, the coefficients as, as and aio could in general be a function of F and III. We propose 
to expand these coefficients in a power series of F and determine them by using DNS data of 
turbulent shear flows 8 . Figure 3 shows the values of the as, as and aio versus the parameter F 
corresponding to the DNS data The scatter is small that indicates that for the above turbulent 
shear flows the III is small such that the coefficients in Eq. (12) are approximately constant. 
Finally, the coefficients as, a$ and aio can be approximately expressed by a truncated fourth order 
power series of F: 


a 5 = -0.1 — 0.5 F + 1.57 F 2 - 2.22 F 3 + 1.07 F 4 
a 8 = -1.6 F + 4.69 F 2 - 5.92 F 3 + 2.62 F 4 

aio = (0.2 + 2.59 F - 7.4 F 2 4- 8.86 F 3 - 3.73 F 4 ) 

Using Eq. (15) the rapid pressure-strain correlation model becomes 


n! 


(i) 

_ 


2q 


— — — (ptkSjk lyfitjbklSfcl) 


2 2 22 2 
+ a 8 (2 b 2 kl S k iS t] + - b 2 kt S k „ + -b 2 kj S kt - jb kt b h S kl + -II S tJ ) 

— -aio (tfiSji + b 2 iS t i - 2b k: bi t S k i — 3b t} b k iS k i) 


+ g (2 + 7as)(6,fcf2jfc -I- b jk Q tk ) — (g a s + -aio)(b 2 iQ^i + b 2 iCl t i) 


(18) 


Figures 4, 5, 6 and 7 show the direct comparisons between different rapid models and DNS data. 
As shown in figures, the models based on the truncated F power series are very encouraging. 
However, it should be pointed out that the model (18) with Eq. (17) is good for the cases where 
\III\ is small. For the cases where \III\ is not negligible, for example, a turbulence undergoes a 
rapid axisymmetric expansion, the coefficients in Eq. (17) must depend on the third invariant III. 
This explains the poor comparison between the model and DNS data 10 for axisymmetric strain 
and plain strain turbulence shown in Figures 8, 9 and 10. 


Improvement of LRR model 

The Launder, Reece and Rodi’s model, Eq. (8), is tensorially linear in 6^ which satisfies the 
basic conditions in Eq. (5), but does not satisfy the realizability condition in Eq. (9). In this 
model, there is only one undetermined coefficient C2 which is shown to be not a constant. We 
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may consider it to be a function of F and then use DNS data to determine its function form. For 
turbulent shear flows, the most important component is the 1-2 component. Therefore, we use 
DNS data of II 12 to deduce the coefficient C 2 . A simple quadratic function is found as follows 

C 2 = 0.38 -0.6 f + 0.5 F 2 

The LRR model with the new coefficient C 2 leads to excellent agreement with all DNS data for 
the 1-2 component (see Figures 4, 5, 6, 7). However, other components are unable to be predicted 
very well. For comparison, another tensorially linear model in kj proposed by Speziale, Sarkar 
and Gatski 11 (SSG) is also included in Figures 4-10. The behavior of the rapid part of SSG model 
is similar to LRR model. 
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Figure 1. The model constant C2 of LRR model 
deduced from the DNS of turbulent shear flow C128W. 







Figure 3. The model coefficients 05, as and aio in the present model 
deduced from the DNS data of turbulent shear flows. 
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Figure 9. Direct comparison betwenn the rapid models and 
the DNS data of the axisymmetric expansion flow EXO. 
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